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The past and recent data analyses gave some hints ol steps in dark energy. Considering the 
dark energy as a dynamical scalar field, we investigate several models with various steps: a step 
in the scalar potential, a step in the kinetic term, a step in the energy density and a step in the 
equation-of-state parameter w. These toy models provide a workable mechanism to generate steps 
and features of dark energy. Remarkably, a single real scalar can cross w = — 1 dynamically with a 
step in the kinetic term. 



PACS numbers: 95.36.+X, 98.80.-k 



I. INTRODUCTION 

Dark energy is a popular explanation for the recent 
acceleration of our Universe. In order to draw a por- 
trait of the dark energy, it is necessary to parameterize 
it and constrain the parameters from observational data. 
By the past and recent data analyses, it was indicated 
that the dark energy could have various steps in its pa- 
rameters. Several years ago, a late-time transition in the 
equation of state (EOS) was studied by Bassett et al 
Very recently, Huang et al [|J reported that the dark en- 
ergy might spring out at a low redshift. However, there 
were few theoretical studies to explain such steps in the 
EOS or in the density of dark energy. 

Our purpose here is to make a step toward this di- 
rection. We will treat the dark energy as a dynamical 
scalar field, and then study a step in the scalar poten- 
tial, in the kinetic term, in the energy density and in the 
EOS respectively. Lacking of enough observational data, 
hitherto the portrait of dark energy is still vague and the 
hints of steps are weak, so we do not attempt to build 
a realistic model in this paper. Instead, to give a vivid 
picture of the mechanism, we will play with simplified 
toy models and choose exaggerated model parameters. 
Because these models are difficult to solve analytically, 
for the most time we will rely on numerical algorithms. 
The models and algorithms can be easily extended and 
refined to give a more realistic description of the dark 
energy. 

Although this could be the first time to systematically 
study cosmic steps in dark energy models with explicit 
Lagrangian, the stepped model is not novel in cosmology. 
It is an old story: possibly steps of the inflaton poten- 
tial have left some fingerprints at the born time of our 
Universe Since the dark energy is more elusive, the 



physical origin of the stepped dark energy field is less 
clear than the stepped inflaton field. 



II. METHODOLOGY 

Before going to specific models, we will describe the 
general framework and numerical methods in some detail. 
Impatient readers can skip directly to the next section to 
get our models (fT0| , (|13p , (|20|) and main results depicted 
in figures. 

In the absence of spatial fluctuations, the Lagrangian 
density of a scalar field minimally coupled to gravity has 
the form 
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where a is the scale factor. Here the function f(<f>) in the 
kinetic term is new. It is positive for quintessence [i[ and 
negative for phantom Later on we will also discuss 
a new model in which /(</>) evolves dynamically from +1 
to —1. In that situation, the single real scalar plays the 
role of quintom 

In a flat universe dominated by dark energy together 
with cold dark matter, if we ignore the contribution of or- 
dinary matter for simplicity, then the evolution dynamics 
is governed by the following Friedmann equations 
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Here p m is the energy density of dark matter, taking value 
p m o at the present time with the scale factor a = ap. The 
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dark energy has an energy density p and a pressure p, 
whose subscripts have been left out for briefness. One 
should not mistake p and p as the total energy density 
and the total pressure. The EOS parameter is defined by 
w = p/p as usual. 

In the above we have used a dot to denote the deriva- 
tive with respect to comoving (physical) time t. For in- 
stance, we have taken the convention of notation = 
d(f>/dt and defined H = a/a. It will be convenient 
to employ the notations x = ln(a/ao) = — ln(l + z) 
and 0' = d(j)/dx, then they give us a useful relation 
= Hcj)' = -H(l + z)(/) tZ . Utilizing 0, we find the 
kinetic energy is 
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The equation of motion 

for the scalar field can be written as 

4*G N f(l + z) [p m0 {l + zf + 4V] 0, 2 
= 8ttGW/(1 + zf [ Pm0 (l + zf + V] <l> tZZ 
+16n 2 G 2 N f 2 (l + zf [p m0 (l + zf + 2V] . 
+4 7 rG Ar /, (l + zf [p m0 (l + zf + V] 4> 2 Z 
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Given explicit form of /(0) and V(<j>), equation © alone 
dictates the evolution of dark energy. But this is a highly 
nonlinear second order differential equation. In most 
cases we have to resort to numerical methods. We will 
take two slightly different schemes, dubbed the linear it- 
eration method and the cubic iteration method, since 
they reform equation ([6]) into linear algebraic equations 
or cubic algebraic equations respectively. Let us elabo- 
rate on them now. 

For numerically evolving equation ^ from an initial 
point Zi to the final point Zf, we partition the interval 
[zi, Zf] into N equal subintervals of width h = Zj — Zj^i = 
(zf — Zi)/N. It is convenient to notate <pj = cf>(zj) for 
short, where j — 0, 1, 2, N. For numerical computa- 
tion, the derivatives can be approximated by finite dif- 
ferences 
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Making use of approximations ([7]) and ©, one can re- 
cast equation ([6]) into a linear algebraic equation with 
respect to <pj+2- Starting with the initial values of cj)j at 
j = 0, 1, 2, 3, we can get the values of all <f>j with j > 3 
iteratively. The linear iteration method gives a unique 



path of <f> for numerical evolutions. Its disadvantage is 
the excessive number of initial conditions. Remember 
that for a second order differential equation we usually 
impose two initial conditions. In practical operation, we 
treat with potentials flat at the initial point 4n, and set 

00 = 01 = 4>2 = ^3 = <j>i- 

Rather than ([5]) one may tend to estimate the second 
order derivative with 
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Now we need only two initial conditions, for instance, 
0o =0i = 0i- This is the starting point of the cubic 
iteration scheme. In such a scheme, instead of a linear 
equation of 0j+2, one has to solve a cubic equation with 
respect to 0j+i- The formula of roots of a cubic equation 
is well known. For our purpose, it is enough to treat 
them in two general categories: (i) if the equation has 
one real root and a pair of imaginary roots, then 0j+i 
takes a value of the real root; («) if all roots are real, we 
determine the value of 4>j+i by minimizing \4>j+i — 4>j\- 
The trick enables us to pick out the smoothest evolution 
path of 0, and fortunately this path is unique in our 
simulation. Even if the path is not unique after applying 
the above trick, one can still find the smoothest path by 
further minimizing |0j+i + <f>j—i — 20j-|, etc. 

The above two methods are operated independently 
in our algorithms. As double checks, they agree with 
each other very well. In fact, the resulting graphs look 
like duplicates. During the numerical simulation, we 
define the reduced Planck mass M p i = \/\/EttGn, the 
fractional density of dark energy f^E = Pde/{3M 2 H 2 ) 
and the present fractional density of dark matter TO o = 
Pmo/{SM 2 H 2 ), where H is Hubble parameter at the 
present time. Moreover, we set f2 m0 = 0.3 and work 
in the unit M p i = 1. In a flat universe, ignoring the con- 
tribution of ordinary matter, one has JlDE + ^m = 1- For 
every specific model below, we evolve equation © from 
Zi = 20 to Zf — with N = 2000 subintervals. Fixing the 
initial conditions and other parameters, if we increase ei- 
ther N or Zi, the change in results is unobservable. This 
confirms the reliability of our methods above and results 
below. 



III. MODELS WITH STEPS 
A. Steps in the Scalar Potential 

So far we have not specified our models. It will be 
done in this section. Firstly, suppose the potential of 
dark energy is broadly flat but has a sudden transition 
near = 6. Then it is natural to model such a stepped 
potential in the following form 
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Figure 1: (color online). Evolution curves of the quintessence 
field with potential (|lUp and / = +1. There are bumps in 
the curves of fractional density and EOS parameter, although 
the bump in fractional density is almost unnoticeable. The 
amplitude and location of bumps can be tuned by changing 
parameters in our model. We set SttGn = 1 in all figures. 



Figure 3: (color online). For stepped phantom with poten- 
tial (|10l) and / = — 1. The dark energy density climbs up a 
step, while EOS parameter has a dip nearby. Similarly to the 
stepped quintessence, the amplitude and location of the steps 
and dips are tunable as we change the model parameters. 
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Figure 2: (color online) . The rise and fall of the kinetic energy 
and potential of a stepped quintessence with potential (|10p 
and / = +1. Both the kinetic energy and the potential are 
normalized by 3Hq, where Ho is the Hubble parameter at 
z = 0. 



For relatively small d, the hyperbolic tangent function 
is a good smooth analytic approximation to a step func- 
tion. The parameter b determines the location of step in 
field space, while c and d control the height and width 
of step respectively. Based on this model, more compli- 
cated ones can be obtained by superposing the steps or 
replacing the constant Vb with other functions. Note here 
c is a dimensionless parameter, but b and d are in unit of 
reduced Planck mass M p i — 1/^/8ttGn- However, in our 
simulation and figures, we will simply set M p i = 1. 
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Figure 4: (color online) . The evolution of the scalar field, dark 
energy density and EOS parameter for the stepped quintom 
model l)13p. In contrast with previous figures, a small change 
in parameter d does not modify the evolution curves signifi- 
cantly. 



We study this potential in the context of quintessence 
with / = +1 and phantom with / = — 1 respectively. 
Since ^deo + ^mo = 1 in a flat universe with negligi- 
ble ordinary matter, it is easy to prove Vq — 3Hy(l — 
^mo)/(l i l c l) f° r quintessence (lower sign) and phan- 
tom (upper sign). The simulation results for stepped 
quintessence are shown in figures [1] and [2] The results 
for phantom are depicted in figure [3] In plotting them 
we have chosen the step location 6 = 0, and the initial 
value of (j> is chosen as fa = —Ad at redshift Zi = 20. 
From the figures we can see a step in the energy density. 
Quintessence falls down the step but phantom climbs up. 
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Although the location of step in field space is fixed by 
b = 0, its location in redshift space is dependent of c and 
d. At the same location there is a feature (a bump for 
quintessence or a dip for phantom) in curves of kinetic 
energy and EOS parameter. Amplitudes of steps and fea- 
tures are determined mainly by parameter c. The initial 
value fa of the scalar field resides well in the flat region of 
the potential because tanh(— 4) ~ —0.999. If we increase 
the absolute value of fa, the steps and features will shift 
to lower redshift region. 

From figures [5] and O one can roughly read the loca- 
tion zt and width At of the step as well as the initial 
density pi and final density pf of the dark energy. In 
data analysis, it is useful to capture these properties by 
parameterizing the dark energy density as 
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This parametrization is applicable when the dark energy 
changes from pi to pf near the redshift zt- Assuming the 
dark energy is decoupled with other ingredients, we have 
the continuity equation p + 3H(p+p) = or equivalently 
(hip)' = —3(1 + w). Corresponding to (JTTJ) , it gives 



(l + z)( Pl - P/ )exp(^f) 
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One may check that the function gives rise to a bump if 
Pi > pf or a dip when pi < pf, in agreement with figures 
|T]and|l 



B. Steps in the Kinetic Term 

Having discussed the steps in the dark energy poten- 
tial, we discuss what will happen if there is a step in the 
coefficient f{fa) of the kinetic term. It is interesting to 
study a new model 



f{fa) — ctanh 
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1 + tanh [ 

1 



(13) 



In favor of the fact that Odeo + ^mo = 1, we can set 
V = 3H§(1 - n m o)/2. In the potential V(fa) there is 
a deep dip, whose depth has been fixed. The coefficient 
/(</>) is non-trivial. It has a step located at <fi = b with 
amplitude c and width d. Again c is a dimensionless 
parameter, while b and d are in unit of reduced Planck 
mass M p i. As mentioned before, we set M p i = 1 in our 
simulation and figures. 

The idea is to settle the potential bottom <f> = in the 
region / > by ensuring ctanh(— b/d) > 0. Given an ap- 
propriate initial condition, the scalar first rolls down to 



the bottom of potential and then moves up, getting closer 
and closer to the point <p — b. However, this model still 
requires some fine-tuning, because it is not always possi- 
ble for the scalar field to pass the sign-inversion point of 
/(</>) after leaving the bottom of potential. We fine-tune 
the parameters, and then numerically evolve the equation 
of motion ([5]). According to the simulation results in fig- 
ure |H both the density and the EOS parameter of dark 
energy decrease near the sign- inversion point. In par- 
ticular, the EOS parameter jumps down abruptly from 
w > — 1 to id < - 1. This is a simple way to cross 
the EOS barrier w = — 1. To make it we need only one 
real scalar field, completely relying on its own dynamics. 
However, one should not regard ([13"]) as more than a toy 
model. There is an infamous instability in any model 
with a wrong-signed kinetic term, e.g., in the phantom 
model Q. We suspect the stepped quintom model here 
suffers from the same problem. Various issues on this 
class of model were explored in HUES]. 



C. Steps in the EOS Parameter 

If there is a step in the EOS parameter, one usually 
parameterizes w in a form [l[ similar to (|lip . However, 
for building a model later, we parameterize it in an al- 
ternative way, 
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where A > 1. With this form, the continuity equation 
can be integrated out as 
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Note there is a duality Wi ^ Wf, \ ^ — A in (fT4|) . We 
are most interested in the following special cases: 



• Wi = — 1; Now the EOS takes the form 
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Wf = —1. In this case, we reexpress the EOS as 

A/[3(l+«0] 
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The equations of state (fT6"]) and (TT7| can be casted into 
a unified form 
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with = a^/67rG w /(l + w c ) and V c = [(1 + 

)/A] 1 / Q . For various choices of parameters, the shape 
of this potential is illustrated in figure [5] The EOS pa- 
rameter w — > w c asymptotically as <f> — > ±oo. Near 
the bottom or top of the potential, it approaches — 1. 
The field rolls down the potential as a quintessence if 
w c > — 1, but rolls up as a phantom when w c < — 1. 
We impose |3cv(l + w c )\ ^> 1 by hand to accomplish the 
sudden EOS transition ([Ml). 



Figure 5: (color online). A rough illustration of potential 
(|20|) with different parameter choices. The arrows indicate 
evolving directions of the dark energy field. For example, 
when a < and w c < — 1, the potential is illustrated by the 
thin solid blue line, and the scalar field rolls up the potential 
from the bottom <f) — to <f> —* ±oo, while the EOS parameter 
W evolves from —1 to w c asymptotically. The dashed blue 
line corresponds to the potential with a < and w c > — 1, 
in which case the scalar field rolls down the potential from 
<f> — » Too to one of the local minimum point, and the EOS 
parameter w evolves form w c to — 1. 



from which one can integrate the continuity equation to 



give 
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It is time for us to reconstruct the potential of ([T]) 
from (|14p in a universe dominated by dark energy. To do 
this, we assume f(<j>) is a nonzero constant and p m = 0, 
W c =/= — 1. For details on reconstructing the potential of 
scalar dark energy, please refer to [DIE* El 13- Making 
use of equations (J5J), ©, (fT%]) and (fT9)) . we finally obtain 



IV. DISCUSSION 



Treating dark energy as a scalar field, we explicitly 
modeled steps in its potential, kinetic term, density and 
EOS, and thus provided a workable mechanism to explain 
the previously and recently claimed dark energy transi- 
tions. To arrive at a realistic model, more experimental 
and theoretical efforts are needed in the future. When 
this work was near completion, a related paper 0] ap- 
peared, where they arrived at models with an implicit 
potential similar to the shape of V((j>) in (fT3|) . As a par- 
tial list, some previouspapers relevant to our work are 

[H d Hzl, EE mm mi- 
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